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Abstract 

We present Monte Carlo data showing the comparison between the parton shower gener- 
ated by the standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-Callan-Symanzik (DGLAP- 
CS) kernels and that generated with the new IR-improved DGLAP-CS kernels recently 
developed by one of us. We do this in the context of HERWIG6.5 by implementing the 
new kernels therein to generate a new MC, HERWIRI1.0, for hadron-hadron interactions 
at high energies. We discuss possible phenomenological implications for precision LHC 
theory. We also present comparisons with FNAL data. 



t Work partly supported by US DOE grants DE-FG02-05ER41399. DE-FG02-09ER41600 and by 
NATO Grant PST.CLG.980342. 



1 Introduction 



With the advent of the LHC, we enter the era of precision QCD, which is characterized 
by predictions for QCD processes at the total precision [TJ tag of 1% or better. At such 
a precision as we have as our goal, issues such as the role of QED [2,3J are an integral part 
of the discussion and we deal with this by the simultaneous resummation of QED and 
QCD large infrared(IR) effects, QED £g> QCD resummation [I] in the presence of parton 
showers, to be realized on an event-by-event basis by MC methods. We stress that, as 
shown in Refs. j3], no precision prediction for a hard LHC process at the 1% level can be 
complete without taking the large EW corrections into account. 

In proceeding with our discussion, we first review our approach to resummation and its 
relationship to those in Refs. [HIE] ; this review is followed by a summary of the attendant 
new IR-improved [7JIE] DGLAP-CS theory [9l[T0] with some discussion of its implications. 
We then present the implementation of the new IR-improved kernels in the framework of 
HERWIG6.5 [11] to arrive at the new, IR-improved parton shower MC HERWIRI1.0. We 
illustrate the effects of the IR- improvement first with the generic 2— >2 processes at LHC 
energies and then with the specific single Z production process at LHC energies. The 
IR-improved showers are generally softer as expected and we discuss possible implications 
for precision LHC physics. We compare with recent data from FNAL to make direct 
contact with observation. Section [5] contains our summary remarks. 

To put the discussion in the proper perspective, we note that the authors in Ref. [T2| 
[T3] have argued that the current state-of-the-art theoretical precision tag on single Z 
production at the LHC is (4.1 ±0.3)% = (1.51 ±0.75)% (QCD) © 3.79% (PDF) © 0.38 ± 
0.26% (EW) and that the analogous estimate for single W production is ~ 5.7%. These 
estimates, which can be considered as lower bounds, show how much work is still needed 
to achieve the desired 1.0% total precision tag on these two processes, for example. This 
point cannot be over-emphasized. 



2 QED(g)QCD Resummation 

In Refs. [UEJE], we have derived the following expression for the hard cross sections in 
the SM SU 2L X U x X SU% EW-QCD theory 



j a = e suM IR (QCED) 1 f d 3 P2 d 3 q 2 tt d?k h -i-r d 3 k' 

exp ^ n\m\ J pS On kj, 
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: By total precision of a theoretical prediction, we mean the technical and physical precisions combined 
in quadrature or otherwise, as appropriate. 
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where the new YFS-style [2] residuals (3 n>m {ki, . . . , k n ; k[, . . . , k' m ) have n hard gluons 
and m hard photons and we show the final state with two hard final partons with mo- 
menta P2, q2 specified for a generic 2f final state for definiteness. The infrared func- 
tions SUMtr(QCED), -Dqced are defined in Refs. [HOE]- This is the simultaneous 
resummation of QED and QCD large IR effects. Eq. ([T|) is an exact implementation of 
amplitude-based resummation of the latter effects valid to all orders in a and in a s . 

Our approach to QCD resummation is fully consistent with that of Refs. [5JE] as 
follows. First, Ref. [15] has shown that the latter two approaches are equivalent. We 
show in Refs. [71 [8] that our approach is consistent with that of Refs. [5] by exhibiting 
the transformation prescription from the resummation formula for the theory in Refs. [5] 
for the generic 2 — > n parton process as given in Ref. [H] to our theory as given for 
QCD by restricting Eq.([T]) to its QCD component, where a key point is to use the color- 
spin density matrix formulation of our residuals to capture the respective full quantum 
mechanical color-spin correlations in the results in Ref. [16] - see Refs. [7,8J for details. 

We show in Refs. [7JE] that the result Eq.flT]) allows us to improve in the IR regime 
the kernels in DGLAP-CS [9j[l0] theory as follows, using a standard notation: 



P%>{z) = C F F YFS { lq )e 



2 "q 



1 + z 
1 - z 



2 



(i - - fMS{i - z) 



Pgg{z) = 2C G F yFS { lG )e¥o{l^ z io + _L_ (1 _ z ya 

+ ^~ 1+7G (1 -z) + z(l - z) 1+ ^) - fabaW - z)}, 
P^{z) = F YFS { lG )e^ l -{z 2 {\ - z)t° + (1 - zfz^}, (2) 
where the superscript "exp" indicates that the kernel has been resummed as predicted by 



2 This should be distinguished from the also important resummation in parton density evolution for 
the "z — > 0" regime, where Regge asymptotics obtain - see for example Ref. p~7l[l8]. This improvement 
must also be taken into account for precision LHC predictions. 
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Eq.fCQ) when it is restricted to QCD alone and where 

7g = Cf— * = -5-, = ^ + — — ( T - ~, 

TV fjQ I IX 6 A 

2 2 1 



fq{lq) = — 



lq lq + 1 7g + 2 

« s 4C G 7g a s C G vr 2 1 

7g = C G —t = —, S G = — + — — (— - -), 

7T Po 2 7T O 2 

jGxiG) qc g F yfs {ig) 70(1 + 7o)(2 + 70) (1 + 7o)(2 + 7g) ' 1 J 

1 1 



2(3 + 7G ) (4 + 7G ) (2 + 7G ) (3 + 7 g) (4 + 7 g) ' 

e -G 7 

iW7) = r(1+7) , C = 0.57721566..., (4) 

where r(w) is Euler's gamma function and C is Euler's constant. We use a one-loop 
formula for a s (Q), so that 

2 

A) = 11 - gn/, 

where n/ is the number of active quark flavors and CV = 4/3 and Cq = 3 are the respective 
quadratic Casimir invariants for the quark and gluon color representations - see Refs. [7J[E] 
for the corresponding details. The results in Eq.([2]) have now been implemented by MC 
methods, as we exhibit in the following sections. 



3 Illustrative Results/Implications 

Firstly, we note that the connection to the higher order kernels in Refs. [19] has been 
made in Ref. [7]. This opens the way for the systematic improvement of the results 
presented herein. Secondly, in the NS case, we find (7j that the n = 2 moment is modified 
by ~ 5% when evolved with Eq.fl2]) from 2GeV to lOOGeV with n f = 5 and A QCD = 
0.2GeV, for illustration. This effect is thus relevant to the expected precision of the 
HERA final data analysis [2U]. Thirdly, we have been able to use Eq.((T]) to resolve 
the violation J2TJ [22] of Bloch-Nordsieck cancellation in ISR(initial state radiation) at 
0(a 2 s ) for massive quarks [23]. This opens the way to include realistic quark masses as we 
introduce the higher order EW corrections in the presence of higher order QCD corrections 
- note that the radiation probability in QED at the hard scale Q involves the logarithm 
ln(Q 2 /m 2 ), and it will not do to set m q = to analyze these effects in a fully exclusive, 
differential event-by-event calculation of the type that we are constructing. Fourthly, the 
threshold resummation implied by Eq.flTJ for single Z production at LHC shows a 0.3% 
QED effect and agrees with known exact results in QCD - see Refs. [U|221[25]. Fifthly, 
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we have a new scheme [8] for precision LHC theory: in an obvious notation, 

a = = E/^^(*0^(*>)»W), (5) 

where the primed quantities are associated with Eq.([2]) in the standard QCD factorization 
calculus. Sixthly, we have |I] an attendant shower/ME matching scheme, wherein, for 
example, in combining Eq.flQ with HERWIG [II], PYTHIA [26], MC@NLO [27] or new 
shower MC's [28], we may use either p^-matching or shower-subtracted residuals 

{Pn,m(ki, ■ ■ ■ , k n ; k[, . . . , k' m )} to create a paradigm without double counting that can be 
systematically improved order- by order in perturbation theory - see Refs. pE]. 

The stage is set for the full MC implementation of our approach. We turn next to 
the initial stage of this implementation - that of the kernels in Eq.(T2]). 



4 MC Realization of IR-Improved DGLAP-CS 
Theory 

In this section we describe the initial implementation of the new IR-improved kernels 
in the HERWIG6.5 environment, which then results in a new MC, which we denote by 
HERWIRI1.0, which stands for "high energy radiation with IR improvement. "El 

Specifically, our approach can be summarized as follows. We modify the kernels in 
the HERWIG6.5 module HWBRAN and in the attendant related modules [29] with the 
following substitutions: 

DGLAP-CS P AB IR-I DGLAP-CS P%» (6) 

while leaving the hard processes alone for the moment. We have in progress [30] the 
inclusion of YFS synthesized electroweak modules from Refs. [31] for HERWIG6.5, HER- 
WIG++ [32] hard processes. The fundamental issue is that CTEQ [33J and MRST 
(MSTW after 2007) [.34] best parton densities do not include precision electroweak higher 
order corrections and such effects do enter in a 1% precison tag budget for processes such 
as single heavy gauge boson production in the LHC environment, as we have emphasized. 

For definiteness, let us illustrate the implementation by an example [35],[36], which 
for pedagogical reasons we will take as a simple leading log shower component with a 
virtuality evolution variable, with the understanding that in HERWIG6.5 the shower 
development is angle ordered [33] so that the evolution variable is actually ~ E9 where 9 
is the opening angle of the shower as defined in Ref. [35] for a parton initial energy E. In 
this pedagogical example, which we take from Ref. [35], the probability that no branching 

3 We thank M. Seymour and B. Webber for discussion on this point. 
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occurs above virtuality cutoff Qq is A a (Q 2 ,Ql) so that 

-dt 

which implies 



dA a (t,Q%) = -^A(t,Q 2 )J2 I dz Y Pba{z) > 



K{Q , Qt) =exp 



The attendant non-branching probability appearing in the evolution equation is 



A(Q 2 ,t) 



t = k\ the virtuality of gluon a. 



The respective virtuality of parton a is then generated with 

A a (Q 2 ,t)=R, 

where R is a random number uniformly distributed in [0, 1] . With (note /?o 
here, where n c is the number of colors) 



<*s{Q) 



2tt 



&o log© 



we get for example 



dz- 



2n 



P qG {z) 
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2 1 



so that the subsequent integration over dt yields 

" Q2 1 dt 2 



Ql 3 t & ln(^) 

t .o2 



Finally, introducing / into Eq.flB]) yields 



cxp 




' 3b 



If we now let A a (Q 2 1 t) = R, then 



Mi) 



2 
3b 



which implies 



t = A 2 



Q_ 

A 2 



3bn 

2\ 



Recall in HERWIG6.5 [TT] we have 
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-n r n. 



-(llric-10), n ; = 5 



-BETAF 



(14) 



(15) 



(16) 



;i7) 



where in the last line we used the notation in HERWIG6.5. The momentum available 
after a qq split in HERWIG6.5 [TT] is given by 

(„ DBETAF 
QLST \ 
^5SJ • (18) 

in complete agreement with Eq. (ll6j) when we note the identifications t = QQBAR 2 , A = 
QCDL3, Q = QLST. 

The leading log exercise leads to the same algebraic relationship that HERWIG6.5 
has between QQBAR and QLST but we stress that in HERWIG6.5 these quantities are the 
angle-ordered counterparts of the virtualities we used in our example, so that the shower 
is angle-ordered. 



Let us now repeat the above calculation for the IR-Improved kernels in Eq.([2]). We 



have 



YFS 



z l (\-z 



i7g 



[l-zfz 



(19) 
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so that 

" iz ^m P ^ = 4fW7 )*" . (20) 

2 » 6oln(fJ)(7G + l)(7G + 2)( 7c + 3) 

This leads to the following integral over dt 

" Q2 dt 4F YF s{lG)e 5G/2 

Ql ~6 ln(^) ( 7g + 1)(7g + 2) ( 7g + 3) 

4F Y F5(7c)e 7G/4 „ A 8.369604402 \ 

Ei II. 



bo (7g + 1) (7g + 2) ( 7G + 3) ^ ' ftolny,) ^ 

Here we have used 



(21) 



, o= T= + f ^^_l, (22) 



2 vr V 3 2 , 

with Cg = 3 the gluon quadratic Casimir invariant. We finally get the IR- improved 
formula 

A a (Q 2 ,t) = exp[-(F(Q 2 )-F(t))], (23) 

where 

FW) = jWTdf^! Ei ( 1, 8 369604402^ 

W ' i-o (To + 1) (7G + 2) (to + 3) \ ioln^) y 

and Ei is the exponential integral function. In Fig. [1] we show the difference between the 
two results for A a (Q 2 ,t). We see that they agree within a few percent except for the 
softer values of t, as expected. We look forward to determining definitively whether the 
experimental data prefer one over the other. This detailed study will appear elsewhere [3Z] 
but we begin the discussion below with a view on recent FNAL data. Again, we note 
that the comparison in Fig. [1] is carried out at the leading log virtuality level, but the 
subleading effects suppressed in this discussion will not change our general conclusions 
drawn therefrom. 

For further illustration, we note that for the q — > qG branching process in HER- 
WIG6.5 [TT], we have therein the implementation of the usual DGLAP-CS kernel as 
follows: 

WMIN = MIN(ZMIN*(1. -ZMIN) , ZMAX* (1 . -ZMAX) ) 
ETEST = (1. + ZMAX**2) * HWUALF (5-SUD0RD*2 , QN0W*WMIN) 
ZRAT = ZMAX/ZMIN 
30 Zl = ZMIN * ZRAT * * HWRGEN ( ) 
Z2 = 1. - Zl 

PGQW = (1. + Z2*Z2) ( 25 ) 
ZTEST = PGQW * HWUALF (5-SUD0RD*2, QN0W*Z1*Z2) 
IF (ZTEST .LT. ETEST*HWRGEN ( 1 ) ) GOTO 30 

where the branching of q to G at z =Z1 occurs in the interval from ZMIN to ZMAX set by 
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Figure 1: Graph of A a (Q 2 ,t) for the DGLAP-CS and IR-Improved DGLAP-CS kernels 
Eqs. (TH| [23]) . Q 2 is a typical virtuality close to the squared scale of the hard subprocess 
- here we use Q 2 = 25GeV 2 for illustration. 

the inputs to the program and the current value of the virtuality QNOW, HWUALF is the re- 
spective function for a s in the program and HWRGEN(J) are uniformly distributed random 
numbers on the interval from to 1. It is seen that Eq. (T2"5]) is a standard MC realization 
of the unexponentiated DGLAP-CS kernel via 

a,(Qz{l - z))P Gq (z) = a s (Qz(l - z)) 1 + (1 ~ (26) 

z 

where the normalization is set by the usual conservation of probability. To realize this with 
the IR- improved kernel, we make the replacement of the code in Eq. (l25l) with the lines 
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NUMFLAV = 5 

BO = 11. - 2 . /3 . * NUMFLAV 
L = 16./(3.*B0) 

DELTAQ = L/2 + HWUALF(5-SUD0RD*2, QNOW*WMIN) *1 . 184056810 
ETEST = (1. + ZMAX**2) * HWUALF(5-SUD0RD*2 , QNOW*WMIN) 

* EXP(0.5*DELTAQ) * FYFSQ (NUMFLAV- 1) * ZMAX**L 
ZRAT = ZMAX/ZMIN 
30 Zl = ZMIN * ZRAT * * HWRGEN ( ) 

Z2 = 1. - Zl ( 27 ) 
DELTAQ = L/2 + HWUALF(5-SUD0RD*2, QN0W*Z1*Z2) *1 . 184056810 
PGQW = (1. + Z2*Z2) * EXP(0.5*DELTAQ) * FYFSQ (NUMFLAV- 1) 
* Z1**L 

ZTEST = PGQW * HWUALF(5-SUD0RD*2, QN0W*Z1*Z2) 
IF (ZTEST .LT. ETEST*HWRGEN ( 1 ) ) GOTO 30 



so that with the identifications 7 g = L, 5 q = DELTAQ, F Y Fs(l q ) = FYFSQ ( NUMFLAV - 1), 
we see that Eg. (127]) realizes the IR- improved DGLAP-CS kernel Pq?(z) via a s (Qz(l — 
z))P^?(z) with the normalization again set by probability conservation. Continuing in 
this way, we have carried out the corresponding changes for all of the kernels in Eq.Q 
in the HERWIG6.5 environment, with its angle-ordered showers, resulting in the new 
MC, HERWIRI1.0(31), in which the ISR parton showers have IR-improvement as given 
by the kernels in Eq. (J6j). □ We now illustrate some of the results we have obtained in 
comparing ISR showers in HERWIG6.5 and with those in HERWIRI1.031 (see footnote 
H|) at LHC and at FNAL energies, where some comparison with real data is also featured 
at the FNAL energy. Specifically, we compare the ^-distributions, p^-distributions, etc., 
that result from the IR-improved and usual DGLAP-CS showers in what followsjf] 

First, for the generic 2— y2 hard processes at LHC energies (14 TeV) we get the com- 
parison shown Figs. |2J [3] for the respective ISR z-distribution and p\. distribution at 
the parton level. Here, there are no cuts placed on the MC data and we define z as 
z = -Eparton/i^bcam where -Ebcam is the cms beam energy and -Eparton is the respective par- 
ton energy in the cms system. The two quantities z and "p\ for partons are of course not 
directly observable but their distributions show the softening of the IR divergence as we 
expect. 

4 In the original release of the program, we stated that the time-like parton showers had been com- 
pletely IR-improved in a way that suggested the space-like parton showers had not yet been IR-improved 
at all. We subsequently introduced release 1.02 in which the part of the space- like parton showers asso- 
ciated with HERWIG6.5's space-like module HWSGQQ for the space-like branching process G — » qq was 
IR-improved. Recently, the remaining un-IR-improved aspect of the space-like branching process, that in 
HERWIG6.5's space-like module HWSFBR, as also been IR-improved in release 1.031. All of the results 
in this paper were obtained using the latter release. 

5 Similar comparisons for PYTHIA and MC@NLO are in progress and we show some results with 
MCONLO below. 
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Energy fraction of parton shower constituents for QCD 2 -> 2 hard parton scattering. 




Figure 2: The ^-distribution (ISR parton energy fraction) shower comparison in HER- 
WIG6.5. 

Turning next to the similar quantities for the 7r + production in the generic 2— >2 hard 
processes at LHC, we see again in Figs. HJ El that the former spectra are very similar in 
the soft regime while the latter spectra are softer in the IR-improved case. These spectra 
of course would be subject to some "tuning" in a real experiment and we await with 
anticipation the outcome of such an effort in comparison to LHC data. 

We turn next to the luminosity process of single Z production at the LHC, where 
in Figs. [6] and [7] we show respectively the ISR parton energy fraction distribution and 
the Z pt distributionwith cuts on the acceptance as 40GeV < Mz, Pt > 5GeV for 
Z — > up, - all lepton rapidities are included. For the energy fraction distribution we again 
see softer spectra in the IR-improved case whereas for the pt distributions we see very 
similar spectra. We look forward to the confrontation with experiment, where again we 
stress that in a real experiment, a certain amount of "tuning" will affect these results. We 
note for example that the difference between the spectra in Fig. [7J while it is interesting, 
is well within the range that could be tuned away by varying the amount of intrinsic 
transverse momentum of partons in the proton. The question will always be which set of 
distributions gives a better \ 2 P er degree of freedom. 

Finally, we turn to the issue of the IR cut-off in HERWIG6.5. In HERWIG6.5, there 



10 




Figure 3: The p^-distribution (ISR parton) shower comparison in HERWIG6.5. 

are IR cut-off parameters used to separate real and virtual effects and necessitated by the 
+- function representation of the usual DGLAP-CS kernels. In HERWIRI, these parame- 
ters can be taken arbitrarily close to zero, as the IR-improved kernels are integrable [7] l8]R . 
We now illustrate the difference in IR cut-off response by comparing it for HERWIG6.5 
and HERWIRI: we change the default values of the parameters in HERWIG6.5 by factors 
of 0.7 and 1.44 as shown in the Fig. [HJ We see that the harder cut-off reduces the phase 
space only significantly for the IR-improved kernels and that the softer cut-off has also 
a small effect on the usual kernels spectra whereas as expected the IR-improved kernels 
spectra move significantly toward softer values as a convergent integral would lead one to 
expectlll This should lead to a better description of the soft radiation data at LHC. We 

6 We note that in the current version of HERWIRI, the formula for a s (Q) is unchanged from that in 
HERWIG6.5 so that there is still a Landau pole therein and this would prevent our taking the attendant 
IR cut-off parameters arbitrarily close to zero; however, we also note that this Landau pole is spurious 
and a more realistic behavior for a s (Q) as Q — > from either the lattice approach [38] or from other 
approaches such as those in Refs. [311110] could be introduced in the regime where the usual formula for 
a s (Q) fails and this would allow us to approach zero with the IR cut-off parameters. 

7 One must note here that the spectra all stop at approximately the same value zq = .00014 — .0016 
which is above some of the modulated IR-cut-off parameters, as the HERWIG environment has other 
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Energy fraction of re* in HERWIG6.5 for QCD 2^2 hard parton scattering. 




Figure 4: The tt + energy fraction distribution shower comparison in HERWIG6.5. 

await confrontation with experiment accordingly. 

We finish this initial comparison discussion by turning to the data from FNAL on 
the Z pt spectra as reported in Refs. [4TJH2]- We show these results, for 1.96 TeV cms 
energy, in Fig. EO For these DO px data, we see that HERWIRI1.031 gives a better fit 
to the data compared to HERWIG6.510 for low pt, (for pt < 8GeV, the x 2 /d.o.f. are 
~ 2.5 and 3.3 respectively if we add the statistical and systematic errors), showing that 
the IR-improvement makes a better representation of QCD in the soft regime for a given 
fixed order in perturbation theory. We have also added the results of MC@NLO |27^l 
for the two programs and we see that the 0(a s ) correction improves the x 2 /d.o.f for the 
HERWIRI1.031 in both the soft and hard regimes and it improves the HERWIG6.510 
X 2 /d.o.f for pt near 3.75 GeV where the distribution peaks. These results are of course 
still subject to tuning as we indicated above. 

built-in cut-offs that prevent such things as a s argument's becoming too small. What the curves in Fig. [8] 
show then are the relative "relative" probabilities for normalized spectra above Zo- 
8 We thank S. Frixione for helpful discussion on this implementation. 
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Histogram of for te* for QCD 2-> 2 hard parton scattering. 
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Figure 5: The ir + ^-distribution shower comparison in HERWIG6.5. 



5 Conclusions 



In this paper we have introduced the first QCD MC parton showers which do not need 
an IR cut-off to separate soft real and virtual corrections. We have shown that spectra 
at both the parton level and at the hadron level are softer in general. In the important 
process of single Z production, these IR-improved spectra show the expected behavior of 
an integrable distribution. The comparison with the DO pt spectrum in the soft regime 
shows that the IR-improvement does indeed improve the agreement with the data. Of 
course, this just sets the stage for the further implementation of the attendant |J] new 
approach to precision QEDxQCD predictions for LHC physics by the introduction of 
the respective resummed residuals needed to systematically improve the precision tag to 
the 1% regime for such processes as single heavy gauge boson production, for example. 
Already, however, we note that our new IR-improved MC, HERWIRI1.031, available at 
|http: / / thep03.baylor.edu \ is expected to allow for a better \ 2 per degree of freedom in data 
analysis of high energy hadron-hadron scattering for soft radiative effects, thereby enabling 
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Energy fraction distribution of parton shower for single Z production. 




Figure 6: The 2:-distribution(ISR parton energy fraction) shower comparison in HER- 
WIG6.5. 

a more precise comparison between theory and experiment. We have given evidence that 
this is indeed the case. Accordingly, we look forward to the further exploration and 
development of the results presented herein. 
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Generated Z transverse momentum. 
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Figure 7: The Z pr-distribution(ISR parton shower effect) comparison in HERWIG6.5. 
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